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Abstract. In J. Phys. A: Math. Gen. 28, 4305 (1995), K. L. Sebastian gave a path- 
integral computation of the propagator of subdiffusive fractional Brownian motion 
(fBm), i.e. fBm with a Hurst or self-similarity exponent H ^ (0, 1/2). The extension of 
Sebastian's calculation to superdiffusion, H G (1/2, 1], becomes however quite involved 
due to the appearance of additional boundary conditions on fractional derivatives of 
the path. In this paper, we address the construction of the path integral representation 
in a different fashion, which allows to treat both subdiffusion and superdiffusion on an 
equal footing. The derivation of the propagator of fBm for general Hurst exponent is 
then performed in a neat and unified way. 



PACS numbers: 05.40.-a, 02.50.Ey, 05.10.Gg 
1. Introduction 

The Langevin representation for Brownian motion has been known for a long time [HE]- 
In one dimension, it can be formally written as 

x{t) = xo+ [ at'w, (1) 

where ^(t) is a Gaussian, uncorrelated noise function: 

mm = s{t-t'). (2) 

The propagator of Brownian motion, which gives the probability of finding the Brownian 
walker at some location x at time t > after having started at Xq at time 0, is a Gaussian 
whose width grows as t^/^: 

G{x, t\xo, 0) = ^ exp (-^^) • (3) 

However, many interesting processes (see [3] for a review) ranging from hydrology [4] 
to finance [5] or tracer transport by turbulent flows [Hj exhibit anomalous diffusion, 
{{x — XqY) oc t^^, with H 7^ 1/2. In many situations the anomalous behaviour is 
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due to spatio-temporal correlations, which translate into long-term time correlations in 
Lagrangian quantities. A generalization of Brownian motion which was proposed to 
include such correlations is known as fractional Brownian motion (fBm) [7]. Although 
slightly different representations can be found in the literature, the one more convenient 
for our purposes is 



where r(-) is the Gamma function and ^(t) is still Gaussian and uncorrelated li. The 
exponent H corresponds to both the self-similarity exponent of the process [7] and the 
Hurst exponent [10] of the incremental process dx{t) := x{t + dt) — x{t). On the one 
hand, the average motion is invariant under the transformation {x,t) {fi^x,fit). On 
the other hand, the increments dx{t) are uncorrelated for H = 1/2 and anticorrelated 
(correlated) for H G (0,1/2) {H G (1/2,1]). Consequently, (jlj) can be used to model 
diffusive {H = 1/2), subdiffusive {H G (0,1/2)), and superdiffusive {H G (1/2,1]) 
transport processes. 

In a beautiful paper [TT], K. L. Sebastian showed that it is very advantageous to 
rewrite fBm in terms of fractional differential operators. 



where oD^" is the Riemann-Liouville fractional integral of order a [T2l [13] . Indeed, by 
combining path integral methods [H] with fractional calculus, he computed the exact 
propagator of (HI) for if < 1/2. However, the case H > 1/2 is not addressed in [TT] 
because the problem becomes quite more complicated from the technical point of view. 
The difficulties may be explained as follows. Sebastian constructs the path integral 
by inverting ([5]) in terms of Riemann-Liouville fractional differential operators. This is 
easy for H < 1/2 because such inversion is well-defined simply imposing the boundary 
conditions x(0) = xq and x{T) = xt- However, if if > 1/2 additional boundary 
conditions involving fractional derivatives at t = are required. Not only is the physical 
meaning of these 'fractional initial conditions' unclear, but also the computation along 
the lines of [Hj results too cumbersome. 

In this paper, we formulate the path integral for fBm using instead S,{t) as the 
integration variable. This choice avoids the need to invert ([5]), and thus the need for the 
additional fractional initial conditions. In our approach the difficulty resides in how to 
manage the boundary conditions of x{t), which translate into a non-local constraint on 
the admissible realizations of the noise. We will see below that the usage of Lagrange 
multipliers allows to do this in an easy and transparent way, thus providing a general 
path- integral derivation of the well-known propagator of fBm [T5] . 

X In [7, 8 fractional Brownian motion is defined by ([4]) except that the lower limit of the integral is 
taken at — oo and the issue of boundary conditions is not present. Our definition is better suited for the 
class of applications we are interested in (e.g. analysis of series of Lagrangian velocities in turbulent 
plasmas). For such physical problems, the stochastic differential equation is naturally formulated in 
a domain t e (io,oo), with finite ^o- An equation analogous to Q was used in 9 , where fBm and 
fractional time processes are compared. 




(4) 



x{t) =Xq + 



{H+l/2) 



(5) 
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2. Path integral calculation of the propagator 

Given a Langevin equation describing the motion of a particle, the propagator 
G{xti T\xq, 0) is defined as the probability to find the particle aX x = xt at time t = T 
if initially, t = 0, it was located at x = xq- Hence, the propagator can be viewed as 
the motion of the particle averaged over all realizations of the noise compatible with 
the boundary conditions x(0) = Xq, x{T) = xt- In the path integral formalism this is 
formulated as follows. The essential object is the probability measure V{^{t))T>^{t) on 
the space of maps ^{t) : [0,T] R. Since C,(t) is uncorrected in time and distributed 
as a Gaussian for each t, this is naturally defined as 

nmmit) = exp (-^ /^ew'dt) vm- (6) 

We want to compute the propagator of ([5]). Firstly, note that the boundary 
conditions x(0) = Xq, x{T) = Xt translate into the following constraint on ^{t): 

oD~^''+'^'k = XT - xo. (7) 
Therefore, the propagator can be written as the expectation value: 

G{xT,T\xo,0) = js(^oDT^''^'^'^^-{xT-Xo))exp(^-^j\itfdt^Vat), (8) 

where the Dirac delta function ensures that we only integrate over maps ^ (t) satisfying 
©• 

We proceed to compute the path integral on the right-hand side of ([8]). Define the 
action 

s[m] = l£at?dt, (9) 

and consider infinitesimal variations on the space of maps ^{t) satisfying ([7]). We denote 
by ^(t) the map which makes the action stationary under such variations. Observe that, 
in particular, ^(t) verifies ([7]). Now, we perform the following change of variables in (|8]): 

at) = at)+vit). (10) 

The constraint ([7]) implies that 

oD^(^+^/^ = 0. (11) 

Using that the action is quadratic in ^{t) and that C,{t) makes it stationary, 

G{XT, T\xo, 0) = e-^[«Wl j 6 (^oDt^^'^'^'^v) e-'^^^'^^Vr^{t), (12) 

whence we deduce that the propagator is of the form 

G(xT,T|xo,0) = /(T)e-^KW] . (13) 

The function /(T) will be determined at the end of the calculation by imposing the 
normalization condition 

G{xt, T|xo, 0)da;r = 1, VT. (14) 



Path integral formulation of fBm for general Hurst exponent 4 

It only remains to compute the stationary points of 5'[^(t)] subject to the constraint 
([7]). It is difficult to do it directly, for the condition ([7]) is hard to implement. The 
technique of Lagrange multipliers saves the day, though. Our problem is equivalent to 
finding the stationary points of 

S[m, >^'\ = \[ ^(^)'d^ + ^ (0^^^''"''^'^^ - i^T - xo)) , (15) 

under arbitrary infinitesimal variations of ^{t) and the Lagrange multiplier A G M. 
From variations with respect to A, 

5^Sm, A] = 6X (o/^-^^+^/'^e - {XT - xo)) , (16) 

we recover, of course, the constraint ([71): 

^^^(H+l/2)^-^^^_^^_ (17) 

Let us perform now variations in ^(t): 

s^it), A] = r ammt + a ,0-^''^"^^^. (is) 

Jo 

Writing explicitly the fractional integral: 

*5S|{(i),A) = £(mm + ^^^^ fiT - t)«-'iHmAt (w) 

Since 5^5 [f(t). A] must vanish for arbitrary 5^{t), we deduce 

The Lagrange multiplier is determined by applying q-D^''^^^^^'* to Eq. (I2T]) and 
using ( fT71l : 



A 



r(if + 1/2)2 

Immediately, we find 



T 



- xo + ,,,,, / (T - t)^^-Mt = 0. (22) 



A T^-^ 

^^-^°+ r2(g + l/2) 2F = °' (23) 



A = -2i7r2(if + l/2)^:^. (24) 



and 



Inserting this in ( 1211) : 

e(t) = 2ifr(if + l/2)^^^(T - tr-^/\ (25) 
Going back to f[T^ we straightforwardly obtain 

G(xT,T|xo,0) = /(r)exp Li7r2(i7 + 1/2)^=^] . (26) 
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And /(T) is obtained from the normalization property, (1141) : 
Hence, finally, 

G{XT, T|xo, 0) = ^^^^ exp [-HT\H + 1/2)^=^) . (28) 
2.1. The kinetic equation 

Once the propagator is known, we can easily derive the kinetic equation associated to 
fBm. The Fourier transform of G{x,t\xo,0) with respect to x is 

G(k, t) = exp —k^ . (29) 

Differentiating with respect to t: 

^2H-l 

dAk, t) = - 2r2 (g + i/2) ^'^^^' 

Finally, we Fourier invert this expression to get the kinetic equation: 

dtGix, t) = (i/ + i/2) ^'^^"^' (^^^ 

which is a diffusion equation with time-dependent diffusivity, as originally deduced in 
[15] from arguments based on the fluctuation-dissipation theorem. The invariance of 
(|3T|) under the transformation (a;, t) — > {fi^x,fit) is manifest. Obviously, when H = 1/2 
we retrieve the standard diffusion equation associated to Brownian motion: 

dtG{x,t) = ]-dlG{x,t). (32) 



3. Conclusions 



In this paper we have computed by path integral methods the well-known propagator of 
fBm for general Hurst exponent, extending the construction of [11], which was restricted 
to the subdiffusive case. The key technical point in our approach is to use the noise as 
path integral variable instead of the coordinate of the particle trajectory. 
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